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Abstract 

We study theoretically the simultaneous, photo-induced two-particle excitations of strongly correlated 
systems on the basis of the Hubbard model. Under certain conditions specified in this work, the corre- 
sponding transition probability is related to the two-particle spectral function which we calculate using 
three different methods: the dynamical-mean field theory combined with quantum Monte Carlo (DMFT- 
QMC) technique, the first order perturbation theory and the ladder approximations. The results are analyzed 
and compared for systems at the verge of the metal-insulator transitions. The dependencies on the electronic 
correlation strength and on doping are explored. In addition, the account for the orbital degeneracy allows an 
insight into the influence of interband correlations on the two particle excitations. A suitable experimental 
realization is discussed. 

PACS numbers: 71.20.b, 32.80.Rm, 33.55.Ad, 79.20.Kz, 68.49 Jk 
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I. INTRODUCTION 



Correlation among electrons are at the heart of numerous phenomena in condensed matter such 
as the metal to insulator transition, the emergence of magnetic and orbital ordering and high- 
temperature superconductivity 13 . Much of today's understanding of the role of electronic corre- 
lation is based on the analysis of the single particle quantities, e.g. the spectral functions, and how 
these compare with experimental data 4 . Among others, a wide spread experimental technique for 
this purpose is the angle-resolved (single) photoemission spectroscopy, ARPES 4 . On the other 
hand, two-particle quantities are essential for the study of important phenomena such as the opti- 
cal conductivity 5 . Two particle properties may be classified in general into those associated with 
the particle-hole, the hole-hole and the particle-particle channels; different techniques are appro- 
priate to access each of these channels. Probably the most studied one of them is the particle- 
hole channel 6 that governs a number of material properties such as the dielectric and the optical 
response 5 . The particle-particle and the hole-hole channels have been much discussed in connec- 
tion with the Auger electron spectroscopy (AES) 7 ~ 18 and the appearance potential spectroscopy 
(APS)^ 9 -. Early experimental works were focused on simple compounds where the two particle 
spectra are well modeled by a convolution the single particle spectra. For AES/APS from corre- 
lated systems several theoretical works^^^— 2kM have been put forward for the evaluation of the 
two-particle spectral functions, mostly based on the Hubbard model 23 25 . Cini and Sawatzky-^ 
obtained in their pioneering works exact results for a completely filled band within the single-band 
Hubbard model. A number of subsequent studies for arbitrary fillings were conducted, mainly us- 
ing the equation of motion method and the ladder approximation. E.g., in the work of Drchal, 
the equation of motion method was employed to calculate the spectral density of the two-particle 
valence bands 9 based on an approximate single-particle spectral function. Other works^^ 6 - uti- 
lize the ladder approximation but differ in their treatments of the single particle quantities. In the 
work of Treglia et al.—, the one-particle spectrum is calculated by evaluating the second order 
perturbation with respect to the Coulomb interaction and with an additional local approximation 
in order to simplify the calculation. Drchal and Kudrnovsky^ii employed the self consistent 
T-Matrix approximation which is valid at a low electron-occupancy. Seibold et al. 21 proposed 
a new approach based on the time-dependent Gutzwiller approximation (TDGA) 22 to calculate 
the electron-pairing; they compared also their results with those of the bare ladder approximation 
(BLA). These works are mostly discussed in connection with AES and/or APS. Recently, an exper- 
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FIG. 1 . A schematics of the one-photon, two-electron (7,2e) experiment. Upon the absorbtion of a VUV 
photon with an energy o; 7 two electrons are excited into the vacuum and simultaneously detected at the 
energies ei, e 2 and the momenta ki, k 2 . 

imental two-particle technique has been developed in which two (indistinguishable) valence-band 
electrons are emitted and detected with well defined momenta ki and k 2 and specified energies 
6i and e 2 upon the absorption of one single (vacuum ultra violet, VUV) photon 27 30 (the method 
is abbreviated by (7, 2e), i.e. one VUV photon in, two electrons out), as schematically shown 
in FigfT] Excitations by a single electron or positron have also been realized and a variety of 
materials ranging from wide band gap insulators to metals and ferromagnets 31 have been inves- 
tigated. The (7, 2e) technique is the extension of ARPES to two-particles; from a conceptional 
point of view one may then expect to access with (7, 2e) the two-particle spectral properties of 
the valence band, as indeed shown below explicitly. A distinctive feature of (7, 2e) is its vital 
dependence on electronic correlation 32 , i.e. the two electrons cannot be emitted with one single 
photon in absence of electronic correlation. The reason for this is the single particle nature of the 
light-matter interaction in the regime where the experiments are performed. Theoretical studies 
concentrated hitherto on weakly correlated systems such as simple metals^—. Consequently the 
two-particle initial state was modeled by a convolution of two single-particle states with the ap- 
propriate energies. The latter were obtained from conventional band-structure calculations based 
on the density functional theory within the local density approximation. Correlation effects were 
incorporated in the construction of the interacting two-particle states of the emitted photoelectrons 
in the presence of the crystal potential. An exception to this approach is the study of (7, 2e) from 
conventional superconductor where the BCS theory was employed for the initial state 38 . While the 
theory reproduced fairly well the observed experimental trends, the previous theoretical formu- 
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lation will certainly breaks down when dealing with strongly correlated materials such transition 
metal oxides or rare earth compounds with partially filled bands. In particular, features akin to the 
metal-insulator transitions are not captured with previous studies. Experiments for such materials 
are currently in preparations. Hence, it is timely to inspect the potential of (7, 2e) for the study of 
strongly correlated systems. 

In the present work, we will present a general theory for the two-particle photocurrent and in- 
spect the conditions under which the experiment can access information pertinent to the particle- 
particle spectral function in the presence of strong electron correlations. In particular, we will 
inspect the particle-particle excitation in Mott systems at the verge of the metallic-insulating 
transition. For the description of the properties we employ the Hubbard model^— and a non- 
perturbative technique, namely the dynamical mean field theory*^!! (DMFT) in combination with 
quantum Monte Carlo technique (QMC) 41 . For the calculations of the two particle Green function 
we will adopt three different ways: The first one is by calculating the single and the two particle 
spectral functions in the loop of DMFT-QMC self consistently. This ensures the fulfillment of the 
sum rules. In the second and the third approaches we basically follow the methods mentioned 
above for the treatment of AES and APS, i.e. we consider the self convolution (first-order pertur- 
bation) and the ladder approximation (LA). We use however the single particle spectral function 
as obtained from DMFT. 

The paper is structured as follows: In section II we present a general expression for the two- 
particle photocurrent and expose its relation to the two-particle Green function. In section III the 
problem is formulated within the two-band Hubbard model and a discussion is presented on how 
to disentangle matrix elements information from the ground state two-particle spectral density. In 
section IV and V we present and analyze the results for the single and two-band Hubbard model 
and compare the results obtained at various level of approximations. Section VI concludes this 
work. 

II. CORRELATED TWO-PARTICLE PHOTOEMISSION 

The (7, 2e) set up is schematically shown in Fig[T] These experiments are conducted in the 
regime where the radiation field is well described classically and the time-dependent perturbation 
theory in the light-matter interaction and the dipole approximations are well justified (low photon 
density and low photon frequency « 50 eV). An essential point for our study is that the operator 
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D N for the photon-charge coupling is a sum of single particle operators, i.e. D N oc J2i=i -^( r i)-Pi 
(in first quantization) where A is the vector potential and p, ; is the momentum operator of particle 
i. This implies that D N cannot induce direct many-particle processes in the absence of inter- 
particle correlations that help share among the particles the energy transferred by the photon to 
one particle which then results in multiparticle excitations. A mathematical elaboration on this 
point is given in 32 and also confirmed below. To switch to second quantization we write A = 
Smm' (E m '{N v - 2)| A. (pi + p 2 ) \E m (N v - 2)) P 2 . The two-particle photocurrent (J), summed 
over the non-resolved initial and final states n and m is determined according to the formula 33,34 



J = Y E E e-^^ I (E m (N v - 2)\A\E n (N v ))\ 2 5 (E - [E m (N v - 2) - E n (N v )}) 

N v ran 

= ^E E ^ En{Nv) Ml m ,M mm n{E n {N v )\Pl\E m iN v -2)) 

N v mn,m'm" 

5(E-[E m (N v -2)-E n (N v )]). (1) 

Here we introduced the short-hand notation M k \ for the matrix elements. The photon energy is 
denoted by E = hio 1 , and (3 is the inverse temperature. Furthermore, «o = 47r 2 a/w 7 , and a 
is the fine structure constant. P 2 = c a cp stands for the (hole-hole) two-particle operator acting 
on the state with N v particles with the energy E n (N v ). Z is the partition function. Under certain 
conditions specified below (the sudden approximation and for high photoelectron energies), the 
variation of the matrix elements, when we vary co> 7 as to scan the electronic states of the sample, 
is smooth in comparison to the change of the matrix elements of P 2 . Furthermore, the diagonal 
elements of are dominant (see below for a justification), i.e. Mm ~ M. In this situation Eq.® 
simplifies to (p is the density operator). 

J = f EE^^I^-i 2 ^^)!^!^^ - 2 ))( E m(N v - 2)\P 2 \E n (N v )) 

N v ran 



5{E-[E m {N v -2)-E n {N v )}) 



J = E / dte^ E ^\E n {N v )\e im Pte- im P 2 {t = 0)\E n (N v ))e' 



El 



N v n 

T'2 



J = ^Jdt*( P P 2 \t)P 2 (t = 0)) e ™ 



J = / dt < Pl{t)P 2 {t = 0) > e iEt . (2) 
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On the other hand, from the spectral decomposition of the two-particle Greens function one 
infers for the two-particle spectral density P(u) the relation 

P{oj) = ^ EE e ^ (%) l^(^ - 2)|P 2 |£ n (i\y>| 2 (l - e-^S(u -E m - E n ). (3) 

N v mn 

Comparing this equation with Eq.© we conclude that under the assumption M k i m M the photon- 
frequency dependence of the two-particle photocurrent is proportional to the two-particle spectral 
density, i.e. 

We recall that the two particle spectral function obeys the sum rule 

J P(Lu)du = (n t n 4 ) . (5) 
A useful auxiliary quantity is partial double occupancy (up to a frequency H) 

/n 
duP(u). (6) 

IH. THEORETICAL MODEL 

The aim here is to explore the potential of the two-particle photoemission for the study of 
the two-particle correlations in matter. To do so we start from the generic model that accounts 
for electronic correlation effects, namely from the doubly degenerate Hubbard Hamiltonian. In 
standard notation we write 2 ^— 

H = } j t a cl aa c jQa + U 2J n ia ^n ial + U'}] ni la n i2cT ', (7) 

ijaa ia too' 

where t a describes hopping between nearest neighbor sites i, j for the orbitals ae(l,2), U, U' stand 
for the intra- and inter-orbital Coulomb repulsion, respectively. The above Hamiltonian does not 
account for the exchange interaction, pairing and spin flip processes. The Hubbard model even for 
a single band provides an insight into a number of phenomena driven by electronic correlations 
such as the metal-insulator transition which cannot be described usually within a static mean field 
theory or within an effective single particle picture such as the Kohn-Sham method within the den- 
sity functional theory. Within the Hubbard model and for the case of infinite connectivity d — Y 00 
the self energy turns local^^. This fact has lead to the development of a new powerful computa- 
tional scheme for the treatment of electronic correlation, namely the dynamical mean field theory 
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(DMFT). For the practical implementation of DMFT it is essential to map the many -body problem 
onto a single impurity Hamiltonian with an additional self consistency relations 44 . Some of the 
possible applications of DMFT have been discussed in Ref.[40], e.g. the long standing problem 
of the metal insulator transition in the paramagnetic phase is described in a unified manner. From 
a numerical point of view, solving the impurity Hamiltonian is a challenging task in the self con- 
sistency of DMFT. For this purpose, quantum monte carlo (QMC) methods are shown to be an 
effective approach which we will follow in the present work with the aim to calculate the single 
and the two particle Green's functions. We note here that since QMC provides only the data for 
imaginary times or equivalently at certain Matsubara frequencies, we need to perform analytical 
continuation to obtain the zero temperatures dynamical quantities. This we do by means of the 
maximum entropy method that we implemented using the Bryan method. A detailed discussions 
on this topic can be found in Ref. 114711 . 



A. The matrix elements 

Now we have to discuss the validity range of the approximation ([3]) that enabled us to assume 
for the matrix elements Mki ~ M. We consider the experiments in the configuration shown in 
FigHJ The photoelectron momenta k x and k 2 are chosen to be large such that the escape time is 
shorter than the lifetime of the hole states. For the description of the photoemission dynamics we 
concentrate therefore on the degrees of freedom of the photo-emitted electrons (which amounts to 
the sudden approximation). The energy conservation laws reads then (cf. FigQ} 

fk)J 1 — CO = E\ + 62 (8) 

where u is the initial (correlated) two-particle energy. The single particle energies Ej are measured 
with respect to the edge of the valence band (or with respect to the Fermi level \i in the metallic 
case). The matrix elements, e.g. M mm i, reduce in the sudden approximation to two particle 
transition matrix elements M^. The high energy final state (with energies £1,62) we write as a 
direct product of two Bloch states (^k) characterized by the wave vectors k x and k 2 , i.e. 

* kl) k 2 (ri,r 2 ) = ^(rOV^ra). (9) 
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1. Intersite ground state correlation 



Correlation effects enters in the initial two particle states. In absence of spin-dependent scat- 
tering (as is the case here) it is advantageous to couple the spins of the two initial states to singlet 
(zero total spin) and triplet (total spin one) states 45 . In the paramagnetic phase and if the two elec- 
trons are not localized on the same sites (they are mainly around R, ; and R,, with i ^ j) the initial 
state is a statistical mixture of singlet and triplet states. The radial part we write then as 49 



*u)(ri, r 2 ) = [<pi{rx - Rj)y3 2 (r 2 - Rj) ± y?i(r 2 - Rj)y2 2 (ri - Rj)) x(|r 2 - r x + Rj - Rj|) 



The "plus" ("minus" sign) stands for the singlet (triplet) channel. We note that since the transition 
operator D 2 is symmetric with respect to exchange of particles, there is no need to anti-symmetrize 
the final state ©. In Eq. (fT0l) the function ipi (ri — Rj) and ipi (r 2 — Rj) are single particle Wannier 
orbitals localized at the sites Rj and R^, respectively. iVj is the number of sites and x(l r 2 — r i + 
Rj — Rj-|) is a (dynamical) correlation factor which we assumed to be dependent on the relative 
distance between the electrons. The part contains correlation effects due to exchange only. 
Due to the localization of the Wannier states around the ionic sites we expect x( I r 2 — r i + Ri — Rj I ) 
to decay with increasing r 1 / 2 (for i ^ j). Since we are dealing with a lattice periodic problem we 
can express the Wannier functions as the Fourier transform of the Bloch states, i.e. cp(r — Rj) = 
i J2^ BZ ■0q( r ) e ~ iq R * 0--BZ stands for the first Brillouin zone). With this relation and exploiting 
the orthogonality of the Bloch states we obtain upon straightforward calculation the following 
expression for the matrix element 



Bloch state ^ qi , i.e. M qi j ki = (^kjA • pi|^ qi ). In deriving the first term of (fTTT) we assumed 
X(|r 2 - n + Rf - Rj\) to vary smoothly with r 1/2 , i.e. %(|r 2 - ri + Rj - Rj\) m x(|Rj - Rj\) 
for i ^ j. For 3D periodic structure the first two terms of Eq.dTTb vanish (momentum and energy 



^xd^-rt + Ri-Rj 



(10) 




(11) 



In this equation M ( 



is the matrix element for the conventional single photoemission from the 
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conservation laws cannot be satisfied simultaneously). Hence, the transition matrix element is 
determined by the third term of (fTTI) . more precisely by the gradient of the correlation factor \- If 
this gradient is smooth on the scale of the variation of \I/ki.k2 and/or then the matrix element 
vanishes all together since ^/ki,k 2 an d ^ffl are orthogonal. Explicitly we find in this case 

l.BZ 



^ l.BZ 

M if ?z — ^ {exp(-iqi • Ri - iq 2 • Rj)5q 2)k2 5 qijkl ±1^2} 



qiq2 



A • (pi + p 2 )x(|r 2 - ri + Ri - R,|) 



(12) 

r2=0=ri 



From this expression we conclude that the matrix elements diminish for decreasing correlation x, 
in fact for i ^ j this contribution to the pair emission is expected to be marginal due to screening. 



2. On site ground state correlation 

The major contribution to the matrix elements is expected to stem from the onsite emission 
Hi = Hj. Only the singlet state is allowed in the single band Hubbard model. To obtain the 
two-particle wave function we assume in line of the Hubbard model that the two electrons scatter 
via a contact potential of strength U when they are on the same site. The wave function reads then 

^(ri,r 2 ) = [y?i(ri - R;)^ 2 (r 2 - Rj) + y?i(r 2 - Rj^Oi - Rj)] x(|r 2 - n|) 

= *£ 0) x(|r 2 -ri|). (13) 

tyffl describes the on-site two electron states that include exchange correlation only. Using only 
yields zero matrix elements as shown above. To obtain an expression for the correlation factor 
x(|r 2 — ri |) (that tends to 1 for U — > 0) we switch to relative R_ and center of mass coordinates 
R_i_. We find that x(|r 2 — ri|) is determined by the integral (Lippmann-Schwinger) equation (xo 
is determined by asymptotic conditions) x(i?-) = Xo + U J d 3 R'_g r (R_, R'_)5( 3 ^(R'_)x(R / _), 
where g r is the retarded Green's function in the relative coordinate. For (TOT) we find then the 
explicit solution 

* w (r ll p 2 ) = *£°)(r ll r 2 )[l + ^ p (r 1 -r 2 ,0)], (14) 

The key point inferred from this relation is that the two-particle transition amplitude increases as 
U increases (tyffl does not contribute to the matrix elements) and it vanishes for Z7 — >• 0. It should 
be noted here that in general U is a dynamical quantity, as evident from its definitions. 
To summarize this section we can say for fixed momenta k 1; k 2 of the photoelectrons and for a 
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given U, the frequency dependence of the two-particle emission, J(cj), is related to the frequency 
dependence of the spectral function P(u). For a given to, the matrix elements vary with U ; they 
contribute a U 2 dependence to J (to). The additional U dependence of J (to) that stems from the 
spectral function will be inspected below. 



B. Two-particle Green's function 

For our purpose we utilize the general expression for the two particle propagator 

47 ),(QV) (q^ m ) = / (^c k , a v'(r)c q _ k , aiCT (r)4_ P)Qi(J (0)4 )Q , (T ,(0)). (15) 

J is a short-hand notation for — a J" Q dre lulmT and ui m = is the Matsubara frequency, T T 
is an ordering operator for r. The local version of the above two particle Green's function or 
the onsite s-wave electron pair which will be directly calculated in the self consistency loop of 
DMFT-QMC is 

Xpp (r) = <T r At(r)A(0)> where A = c t c ; . (16) 

The evaluation of the two particle propagator may be performed with the aid of the perturbation 
expansion using the standard diagrammatic theory by selecting the diagrams appropriate for the 
physical problem at hand. For the Hubbard model with the short range interaction we utilize the 
ladder-type diagrams. For the single-band Hubbard model (an extension to the multi-orbital case 
is straightforward), the two particle propagator reads 

X PP (q, iu) m ) = -- ^2 £( k ' iu n)G(q - k, iu m - iu n )T(k, q, iu m ). (17) 

kiu„ 

We selected the ladder diagrams and summed to all orders. Since in our model the Coulomb 
interaction is static and independent of the wave vector, the vertex function T reads 

r(k, q, iOJ m ) = 1 - j ^2 £(P> ip 'n)G(Q - P> iuj rn ~ ^n) r (P> ^ iu 'n)i ( 18 ) 

meaning that the right hand side of this relation is independent of k— . Thus we obtain 

w4,u= l-«^) <9) 

where 

X(q, iu m ) = -~n^2 ^^ p ' iv n)Q{^. ~~ P» iu} ™ ~ iv n ) (20) 
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is the two particle Green's function expressed in terms of the full single particle Green's function. 
Performing standard analytical continuation and evaluating the imaginary part of the two particle 
Green's function one arrives at the following expression for the two particle spectral function 



P(co)=lm[ Xpp (q,u + i5)]. (21) 

In order to evaluate the above equation, it is sufficient to calculate the imaginary part of the two 
particle propagator and analytically continue it to real frequencies. This yields 

/oo poo 
dv \ deD(e)[A(e,p)A(-e,u-u)(l-f(p)-f(u-u)}, (22) 
-oo J —oo 

where stands for the imaginary part of xi^), f(u) is the Fermi distribution function and 

A(e,u) = — ^Im U) _ e i E ( aj ) is the full interacting single particle spectral function. D(e) is the 
free density of states and ^ is a constant. The real part of the two particle vertex is obtained 
via the Kramers-Kronig relation, which follows from the causality condition. For the case of the 
degenerate Hubbard model, it is straightforward to extend the above formulation where now each 
Green's function contains the composite orbital spin index, a = (a, a) 

The two particle propagator x a ' a ' (q, iu m ) reads in this case 

1 

1 



\ ( q l^„, ) -- ^ ^(P' iu n)G a \<i - P, iu m ~ iVn)- (24) 



,a,a 



IV. SINGLE BAND HUBBARD MODEL 

An essential ingredient for the calculation of the two particle Green's function (1221) . is the single 
particle spectral function. The results of the spectral function obtained from DMFT-QMC method 
for the single band Hubbard model at half filling and away from filling are presented in Figf2l In 
the QMC calculation, the Hirsch-Fye 41 algorithm is employed with the following parameters: The 
energy scale is set by the bandwidth W = 1. For the temperature we choose T/W = 0.05, the 
increment of time slices is Ar <=0.5. DMFT-QMC calculations are performed for the param- 
agnetic phase and using the Bethe lattice for the free state density in the self consistency loop. 
At half filling (the left panel of Figf2]) the quasi-particle peak at the Fermi energy is the domi- 
nant feature in the single particle spectra in the weak coupling interaction signifying a metallic 
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FIG. 2. The results of the DMFT-QMC for the frequency dependence of the single-particle spectral function 
A(ui) for the single-band Hubbard model at half filling for various interaction strengths U (left panel) and 
for various electron occupancy at U/W = 3 (right panel); for all cases we choose W = 1. 

behavior; the carriers are itinerant and a Fermi liquid picture is appropriate. With an increasing 
strength of electronic correlations, localization sets in accompanied by a gradual disappearance 
of the quasi-particle weight and the formation of a pseudogap. Electron transfer between the two 
bands may occur, albeit its probability is smaller than that in the previous case. As the coupling 
strength further increases, the gap fully develops indicating an insulating state. 

The role of the double occupancy we inspect by studying the quantity (n^ng,) calculated in the 
DMFT-QMC loop. Evolving from the weakly interacting (metallic) case to the strongly interac- 
tion (insulating) phase the double occupancy is reduced 40 , for more energy is required to overcome 
the stronger repulsion whenever forming the double occupation. The influence of dopant concen- 
tration on the MIT is demonstrated by doping the insulating phase as depicted in the right panel 
of Figj2] Contrasting with the results at half filling with an interaction strength U /W = 3, the 
spectral function in this case shows a resonance peak at low energies testifying that the system 
attains again a metallic character. This is because the doping enhances the number of holes which 
in turns increases the itineracy such that the electron can hop from one site to the other. 

Having commented on the generic single particle properties of the single band Hubbard model 
for Mott systems, we turn now to the discussion of the particle-particle spectral function. For small 
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FIG. 3. Two particle spectral function as function of the correlated, two-particle initial energy u) for vari- 
ous interaction strengths. Calculations are performed within the self consistency scheme of DMFT-QMC 
method. The large scale figure is shown in inset (b). Inset (a) shows the integrated spectra using Eq.©. 
Note however, that the energy is measured with respect to the uncorrected two-particle Fermi energy 2/i. 

U/W one obtains an intense peak that lies close to u/W = 0. The origin of such features can be 
inferred from the structure of the single particle spectral function: P in this case is well modeled 
by a convolution of two single-particle spectral functions. A small increase of U/W leads to a 
reduction of the spectral weight which shifts the peak to higher energeis (far from u/W = 0). As 
the interaction strength further increases, the spectral weight decreases significantly signaling a 
reduction of double occupation. This argument is supported by the results of the integrated spectra 
depicted in the inset of Figj3] In addition to the reduction of the spectral weight, one also observes 
the formation of a gap in the low energy regime (near to the zero frequency) for strong interaction. 
This two particle gap resembles the one that appears in the single particle spectra (cf. Fig© which 
is the usual indicator for an insulating state. We argue here that this is also a signal for the system 
in the insulating state from the point of view of particle-particle excitations. As already pointed out 
above, the reappearance of the low energy resonance as a function of the doping is a signal for the 
metallic character and the associated behavior of the single particle spectral function. The same 
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FIG. 4. The same as Figf3]with the same notation, however we inspect here the two particle spectral function 
away from half filling for U /W = 3 and for various electron occupancies. Inset shows the integrated spectra 
according to equation ©. 

pattern is also observed in the two particle spectral function where the strongest peak occurs in the 
lowest electron occupancy and decreases as the Mott insulating phase is approached. Thus the two 
particle spectra also highlight the contribution of holes to the double occupancy probability, and it 
is clearly supported by the results for the integrated spectra (see inset of FigjU). 

A. Relation to the (7, 2e) experiments 

To connect the results of Figj3] to the (7, 2e) signal it is decisive to recall the statements of 
equations (1811 4b : The correlated two-particle initial energy that appears in Eq.® and which is 
scanned in Figj3j is in the uncorrected case merely the sum of two single particle energies tOi 
{Uuncor — oj% + oj 2 ), i-e. in the metallic uncorrelated case we expect some spectral weight around 
to = in Figj3] For a finite U, i.e. for a correlated system one requires more energy to compensate 
for the repulsion of the Coulomb interaction. This is the reason for the shift of the two-particle 
peak in Figj3]with increasing U. The same can be observed in the single particle spectral function 
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where the distance between the two Hubbard band is approximately on the order of U /W. The 
tendency of larger spectral density with decreasing U is not reflected in the (7, 2e) signal J. In 
fact, the opposite will occur. The reason for this is that according to eqs. (PT41I21) J is proportional 
to the product of the matrix elements and the spectral function. On the other hand the matrix 
elements decrease with U (cf. eqs. fll 11141) ). and in fact vanishes for U —s- counteracting against 
the trend with U of the spectral function P (cf. Figj3]). We stress however, that the results shown 
in Figf3]are still relevant to the (7, 2e) measurements in that, for a given U, the matrix elements 
are hardly dependent on oj. 

B. Comparison between different approaches to the two-particle spectral functions 

To inspect the role of the ladder diagram summation (i.e. Eq.(|T9T) with results in FigJ5](b)), we 
compare with the results (shown in Figj5](a)) of the first-order approximation using (|22)) (i.e. with 
the convolution of the single particle spectra). The results of the first order approximation show a 
smooth, broad Gaussian-like feature in the spectra for all interaction strengths. This is due to the 
self-convolution that tends to wash out the character of the original function. The presence of a 
gap in the two-particle spectra highlights the difference between the weak and the strong coupling 
interactions in agreement with the previous result of DMFT-QMC and with the same energetic 
origin as discussed above. That this correct energetic shift is reproduced by this simple scheme 
is the result of using an accurate single particle spectral function. Another point is the evolution 
of the two particle spectra from the weak through the strong coupling limit and the associated 
behaviour of the spectral weight. In the scheme used in Fig|5l the weight seems to be comparable 
for all values of the interaction strengths except for U/W = 2 which originates from the low 
shoulder in the spectra of Figj2] The reduction of the spectral weight is related to the probability 
of the double occupancy. It is then conceivable to infer that this scheme violates the sum rule for 
the two particle spectral function (which is dictated by the double occupancy, see Eq.©. This 
is endorsed by the results for the integrated spectra shown in the inset of FigJ5](a). The shift to 
higher frequencies is due to the presence of the gap. No clear suppression is observed as in Figj3] 
and FigjH 

Having obtained the imaginary part of the first order approximation we inspect the influence of 
the ladder diagram summation on the two particle spectra. The results are presented in FigJ5](b). 
In contrast with previous results obtained in the first order approximation, the spectra delivered 
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FIG. 5. The frequency dependence of the two-particle spectral function at half filling, calculated with the 
first-order perturbation approximation (a) and with the full ladder approximation (b). Various curves cor- 
respond to different interaction strengths. Inset shows the integrated spectra. The single particle quantities 
are obtained from DMFT Same notation and units as Figj3] 



by DMFT-LA are non-uniform with smooth broad feature and a satellite peak. For the weak 
interaction strength, the two particle spectra hardly depends on the Coulomb interaction strength. 
As before no clear reduction of the spectral weight is observed. Interesting features in the DMFT- 
LA scheme emerge at higher interaction strengths, which from the point of view of the single 
particle spectra, is already the regime of the insulating phase. Instead of suppressing the spectral 
weight, the increases of the coupling interaction strength results in a narrow satellite peak. The 
integrated spectra depicted in the inset of Fig. [5](b) shed some light on this result. The integrated 
spectra within the ladder approximation exhibit a suppression of the weight for higher frequencies 
in contrast to results of the first-order approximation. We remark that in the ladder approximation 
the suppression of the integrated spectra is not related to a diminishing of the weight of the two 
particle spectral function but is associated with the width of the spectra that become narrow as 
the interaction increases. All in all we can conclude for these results that the (7, 2e) technique is 
the appropriate tool for testing the validity of approximate schemes for the two-particle Green's 
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FIG. 6. The same as in Fig. ((5]), however we inspect here the role of varying the electron occupancy n 
at an interaction strength of U/W = 3. (a) shows the first-order approximation results whereas in (b) the 
predictions of the ladder approximation are plotted. 

The two-particle spectral function away from half filling is depicted in Figj6]for various occu- 
pancies and for U/W = 3; calculations are performed within the first-order approximation and 
within the ladder approximation. No gap formation in the two particle spectra takes place. This is 
consistent with the behaviour of the single particle spectral function for which the hole doping of 
the insulating phase stimulates the formation of quasi- particles. In the first-order approximation, 
one obtains the usual broad Gaussian-type structure that diminishes as a function of the dopant 
concentration . A somewhat similar situation is also observed for the results of DMFT-LA. In 
the latter, however, one observes an intense low-energy peak in the case close to half filling. The 
peak decreases as the doping increases. The results of both these approaches arein contrast to 
those obtained via DMFT+QMC where the largest spectral weight is obtained for the high doping 
concentration. Therefore, these results do not reflect the fact that an additional doping leads to an 
increase in the double occupancy which is clearly supported by the sum rule results plotted in the 
inset of FiglH Here one observes that the spectral function at the maximum value of the doping 
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obtains the smallest spectral weight. A similar finding has been observed in reference where the 
bare ladder approximation (BLA) has been utilized. In their result, the decrease of the electron 
occupancy also increases the peaks in the spectra, which they assume to be a violation of the two 
particle sum rule. On the other hand, by using the time-dependent Gutzwiller approximation the 
opposite situation occurs: The two particle spectral weight diminishes as the Mott insulating phase 
is approached, which is in line with what we obtained above within the DMFT-QMC. 



V. TWO BAND ISOTROPIC HUBBARD MODEL 



The single band Hubbard model on which we based our above discussion, is useful for systems 
with only a single band being close to Fermi energy. To inspect the role of the orbital degrees 
of freedom, which is known to be important for the properties of strongly correlated systems, a 
multi-orbital model is needed. It is the aim of this section to study the influence of the orbital 
degeneracy on the single and two particle spectra. 

The results for the single particle spectral function within the two band Hubbard model are 
presented in FigjTJ The results are similar to those obtained within the single band Hubbard model 




4 -2 2 4 
coAV 



FIG. 7. The DMFT-QMC results for the frequency dependence of the single-particle spectral function A(uj) 
of the two-band, isotropic Hubbard model at half filling. Various curves corresponds to different interaction 
strengths U in units of the band width W (here W = 1). The insets show the orbitally resolved spectral 
functions for the first (upper inset) and the second band (lower inset). 
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FIG. 8. The two-particle spectral function of the degenerate Hubbard model at half filling as a function of 
the two-particle initial energy u measured in units of W. The calculations are performed by the DMFT- 
QMC method including total bands (a) and inter-band (b). The interaction strengths U are varied. 

(cf. Fig 01)- The metallic phase shows an intense quasi-particle peak that diminishes as the coupling 
interaction becomes stronger. The formation of the gap for a high interaction strength shows the 
existence of the insulating phase in this degenerate system. An essential point that distinguishes the 
Mott transition in the single band from the degenerate band case is the value of the critical coupling 
necessary to obtain a dip in the spectral function. This behavior is well documented in the works 
of reference 48 employing the Gutzwiller approximation. There, a relation has been established 
between the critical coupling and the orbital degeneracy. From the orbitally resolved spectral 
function depicted in the embedded figures, one also learns that each band undergoes the same 
transition from metal into insulating phase. For anisotropic bandwidth each band undergoes an 
independent metal insulator transition, a behaviour coined as the orbital selective Mott transition 46 . 

The results of DMFT-QMC calculations for the two particle spectral function are illustrated in 
Figj8]that contains the two spectral functions of the total band (a) and the interband (b). From 
Figj8]we see that a small increase of the Coulomb interaction in the weak coupling regime hardly 
affects the overall spectral weight. Furthermore, increasing the interaction strength leads however 
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FIG. 9. The same as in Fig.® for the total bands, however we show in (a) the results of the first order 
perturbation approximation and in (b) those of the ladder approximation. 

to the reduction of the spectra as well as to a shift of the dominant peak to higher energy. 

For the case of interband spectra, there is a clear signal of the spectral weight reduction already 
in the metallic case. As the insulating phase is approached, the two particle spectra show a double- 
peak structure. The two particle spectra obtained by means of the first order approximation as well 
as by the ladder approximation are shown in Figj9] (a) and Figj9] (b) respectively. The behavior 
of the two particle spectra in the single band Hubbard model obtained within the same scheme 
(see FigHJ) (e.g. the gap existence, absence of spectral weight reduction) is also observed in the 
present case. In the metallic case however there are new features predicted by both approximations 
namely a double peak structure that disappears in the insulating phase. Other notable features such 
as the increase of the weight as the coupling strength increases are present in the results of both 
methods. The integrated spectra of the degenerate model indicate a violation of the sum rule 
for the two particle spectra by both the first order approximation and the ladder approximation. 
From the three scheme: QMC-DMFT, first order and ladder approximations, the DMFT-QMC 
methods provides the more reasonable predictions which practically always obey the sum rule as a 
constraint on the two particle spectral function. This is because, both the single and the two particle 
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FIG. 10. The same as in Fig.® for the interband Hubbard model at half filling. The figure shows the results 
of the first order approximation (a) and of the ladder approximation (b). 

propagators are calculated on an equal footing in the self consistency DMFT. An accurate single 
particle approach when formulating the two particle propagator 17 , does not however guarantee the 
fulfillment of the sum rules. The use of an accurate approach in the single particle spectra captures 
however pertinent features such as the gap opening in the insulating state which is also observed 
in the result of DMFT-QMC. 



VI. FINAL REMARKS AND CONCLUSIONS 

We shall now comment on the possible implementation of our proposal. It is a widely accepted 
wisdom that the single band Hubbard model can be employed to explain the results of single par- 
ticle or particle-hole properties of vanadium sesquioxide V2O3. It is shown 40 that the variation 
of the Coulomb repulsion U is realized by changing the chemical composition or applying hy- 
drostatic pressure. We therefore also argue in this respect that the two-particle properties as we 
have presented here can be accessed in the similar manner. For the two band Hubbard model, the 
results can be again implemented to describe the physics of V 2 3 . In this case one can investi- 

21 



gate to role of the orbital degrees of freedom. The inclusion of orbital degrees of freedom allows 
the applications of our model to wider class of systems. As for the experimental geometry, our 
theory is limited by the fact that the employed self energy is local. Thus, our predictions are best 
tested by fixing in Fig{T]the momenta of the detected electrons. What should be varied is then the 
photon energy w 7 . Since the momenta and the energies of the detected electrons are fixed by the 
experiment, the only quantity which is scanned is the initial correlated two-particle energy u. A 
simulation of the experiments for varying momenta and fixed u requires an explicitly non-local 
self-energy which goes beyond the validity of the present model. 

To summarize, in this work we explored the potential of two-particle photoemission for the 
study of two-particle correlations in correlated systems. We identified the conditions under which 
the two-particle photocurrent is related to the two particle spectral function. Calculations have 
been performed within the framework of the single and the two band Hubbard model. We per- 
formed calculations and compared the results of three different schemes DMFT-QMC, the first- 
order perturbation and the ladder approximations based on the DMFT single particle spectra. In 
the single band case, the two particle spectral function evaluated with DMFT-QMC is shown to be 
dependent on the double occupancy in the system. As for the single particle spectral function, an 
increase in the electronic correlation strength results in suppression of spectral weight of the two 
particle spectra and in an opening of a gap near zero two-particle frequency. The first-order per- 
turbation and ladder approximation calculations are qualitatively different from the DMFT-QMC 
predictions. A finding that can be directly tested by two-particle photoemission spectroscopy. The 
inclusion of the orbital degeneracy brings about an increase of the critical coupling and additional 
interband contributions to the spectra; these features should also be distinguishable by two-particle 
photoemission experiments. 
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